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Abstract 



Fix k > 1, and let I (I), I > 1, be a sequence of fc-dimensional vectors of multiple Wiener- 
Ito integrals with respect to a general Gaussian process. We establish necessary and sufficient 

■ conditions to have that, as I — > +oo, the law of 1(1) is asymptotically close (for example, 
in the sense of Prokhorov's distance) to the law of a ^-dimensional Gaussian vector having 
the same covariance matrix as 1(1). The main feature of our results is that they require 
minimal assumptions (basically, boundedness of variances) on the asymptotic behaviour of 
the variances and covariances of the elements of I (I). In particular, we will not assume that 
the covariance matrix of I (I) is convergent. This generalizes the results proved in Nualart 
and Peccati (2005), Peccati and Tudor (2005) and Nualart and Ortiz-Latorre (2007). As 

£S) | shown in Marinucci and Peccati (2007b), the criteria established in this paper are crucial 

■ in the study of the high-frequency behaviour of stationary fields defined on homogeneous 
spaces. 

[ — ■ Keywords and Phrases: Central Limit Theorems; Gaussian Approximations; High- 
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■ 1 Introduction 

Let U (I) = (U\ (/),..., Uk (I)), I > 1, be a sequence of centered random observations (not 
necessarily independent) with values in R fc . Suppose that the application I \— > KUi (I) 2 is bounded 
for every i, and also that the sequence of covariances q (i,j) = EC/j (/) Uj (I) does not converge 
as I — > +oo (that is, for some fixed i ^ j, the limit lim^oo q (i,j) does not exist). Then, a 
natural question is the following: is it possible to establish criteria ensuring that, for large I, the 
law of U (I) is close (in the sense of some distance between probability measures) to the law of 
a Gaussian vector N (/) = (JVi (I) ,...,N k (I)) such that EN { (I) Nj (I) = EU { (I) Uj (I) = q (i, j) ? 
Note that the question is not trivial, since the asymptotic irregularity of the covariance matrix 
c\ (•,•) may in general prevent U (I) from converging in law toward a /c-dimensional Gaussian 
distribution. 

In this paper, we shall provide an exhaustive answer to the problem above in the special case 
where the sequence U (I) has the form 



u(0 = i(0 = (i dl (// 1) ),-,^(// fe) )) ; «>i, (i) 
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where the integers di, dt > 1 do not depend on I, I d . indicates a multiple stochastic integral 
of order dj (with respect to some isonormal Gaussian process X over a Hilbert space f) - see 
Section [5] below for definitions), and each 6 i}®^, j = l,...,k, is a symmetric kernel. In 
particular, we shall prove that, whenever the elements of the vectors I (I) have bounded variances 
(and without any further requirements on the covariance matrix of 1(1)), the following three 
conditions are equivalent as I — > +oo: 

(i) 7 (C (I (/)) , C (N (/))) — > 0, where C (•) indicates the law of a given random vector, N (I) is 
a Gaussian vector having the same covariance matrix as I (I), and 7 is some appropriate 
metric on the space of probability measures on W k ; 




(iii) For every j = 1, k and every p = 1, dj — 1, the sequence of contractions (to be formally 
defined in Section EJ) ® p i^, / > 1, is such that 



Some other conditions, involving for instance Malliavin operators, are derived in the sub- 
sequent sections. As discussed in Section [5l our results are motivated by the derivation of 
high-frequency Gaussian approximations of stationary fields defined on homogeneous spaces - a 
problem tackled in [9] and [TO] . 

Note that the results of this paper are a generalization of the following theorem, which 
combines results proved in [13], [T3] and [15]. 



Theorem 0. Suppose that the vector I (I) in (Q]) is such that, as I — > +00, 



EI di (jf) I d . (/<») -» C (i,j) , l<i,j<k, 

where C = {C (i,j)} is some positive definite matrix. Then, the following four conditions are 
equivalent, as I — ► +00: 



1. I (I) A N (0, C), where N (0, C) is a k-dimensional centered Gaussian vector with covari- 
ance matrix C; 

2. Relation {IP takes place for every j = l,...,k and every p = 1, dj — 1; 



3. For every j = 1, ...k, E J d . /, 



r0") 



3C (j,j) 2 ; 



4. For every j = 1, ...A;, 



(j) 



p(j) 



dj in L 2 , where D 



Id, f> 



p(j) 



Malliavin derivative of I / z J (see the next section). 



denotes the 



The equivalence of Points 1.-3. in the case k = 1 has been first proved in p3] by means of 
the Dambis-Dubins-Schwarz (DDS) Theorem (see Ch. V]), whereas the proof in the case 



2 



k > 2 has been achieved (by similar techniques) in [15]; the fact that Point 4. is also necessary 
and sufficient for the CLT at Point 1. has been recently proved in [13], by means of a Malliavin 
calculus approach. For some applications of Theorem (in quite different frameworks), see e.g. 

m, 0, 0, m or mi. 

The techniques we use to achieve our main results are once again the DDS Theorem, com- 
bined with Burkholder-Davis-Gundy inequalities and some results (taken from [4, Section 11.7]) 
concerning 'uniformities' over classes of probability measures. 

The paper is organized as follows. In Section [2] we discuss some preliminary notions con- 
cerning Gaussian fields, multiple integrals and metrics on probabilities. Section [3] contains the 
statements of the main results of the paper. The proof of Theorem [T] (one of the crucial results 
of this note) is achieved in Section 01 Section [5] is devoted to applications. 

2 Preliminairies 

We present a brief review of the main notions and results that are needed in the subsequent 
sections. The reader is referred to [6] or |12|. Ch. 1] for any unexplained definition. 

Hilbert spaces. In what follows, the symbol indicates a real separable Hilbert space, with 
inner product (•, •}« and norm ||-||^. For every d > 2, we denote by fj® 2 and fj 02 , respec- 
tively, the nth tensor product of fj and the nth symmetric tensor product of $). We also 
write Sj m = Sj &1 = Sj. 

Isonormal Gaussian processes. We write X = {X (h) : h G fj} to indicate an isonormal 
Gaussian process over fj. This means that X is a collection of real- valued, centered and 
(jointly) Gaussian random variables indexed by the elements of fj, defined on a probability 
space (p., J 7 , P) and such that, for every h, h' G fj, 

E[X(h)X(h')] = <M%- 

We denote by I? (X) the (Hilbert) space of the real-valued and square-integrable func- 
tionals of X. 

Isometry, chaoses and multiple integrals. For every d > 1 we will denote by 1^ the isom- 
etry between $) &d equipped with the norm V"! and the dth Wiener chaos of X. 

In the particular case where f) = I? (A,A,fi), (A, A) is a measurable space, and fj, is a 
a-finite and non-atomic measure, then fi &d = L 2 S ( y A d , A® d , is the space of symmet- 
ric and square integrable functions on A d and for every / G Sj &d , 1^ (/) is the multiple 
Wiener-Ito integral (of order d) of / with respect to X, as defined e.g. in [121 Ch. 1]. It 
is well-known that a random variable of the type Id(f), where d > 2 and / ^ 0, cannot 
be Gaussian. Moreover, every F G L 2 (X) admits a unique Wiener chaotic decomposition 
of the type F = E (F) + X^2=i (/d)> where fd G $) &d , d>l, and the convergence of the 
series is in 1? (X). 

Malliavin derivatives. We will use Malliavin derivatives in Section [3l where we generalize 
some of the results proved in [13]. The class S of smooth random variables is defined as 
the collection of all functionals of the type 

F = f(X(h l ),...,X(h m )), (3) 
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where hi,...,h m G f) and / is bounded and has bounded derivatives of all order. The 
operator D, called the Malliavin derivative operator, is defined on S by the relation 

M a 
O 



DF = Y J Q^J{hu-,h m )h i , 

where F has the form ([3|). Note that DF is an element of I? (O; S)). As usual, we define the 
domain of D, noted D • , to be the closure of S with respect to the norm H-F^ 2 — E (-F 2 ) 
+E \\DF\\^. When F G D 1 ' 2 , we may sometimes write DF = D [F] , depending on the 
notational convenience. Note that any finite sum of multiple Wiener-Ito integrals is an 
element of B ,2 . 

Contractions. Let {e^ : k > 1} be a complete orthonormal system of S). For any fixed / G Sj &n , 
g G S) Qm and p G {0, n A m}, we define the pth. contraction of / and g to be the element 
of ^ n+m ~ 2p given by 

oo 

i\,...,ip=\ 

We stress that / (& p g need not be an element of ^0 n + m - 2 P. \y e denote by f® p g the 
symmetrization of / ® p g. Note that / (g>o 9 is just the tensor product f ® g oi f and g. If 
n = m, then / ® n g= {f,g)^ n . 

Metrics on probabilities. For k > 1 we define P (M fc ) to be the class of all probability mea- 
sures on IR fc . Given a metric 7 (•, •) on P (M fc ) , we say that 7 metrizes the weak convergence 
on P (M fc ) whenever the following double implication holds for every Q G P and every 
{Qz : I > 1} C P (M fe ) (as I — > +00): j(Qi,Q) — ► if, and only if, converges weakly 
to Q. Some examples of metrizing 7 are the Prokhorov metric (usually noted /?) or the 
Fortet-Mounier metric (usually noted f3). Recall that 

p (P, Q) = inf{e > : P (A) < Q(A € ) + e, for every Borel set A C (4) 

where A e = {x : \\x — y\\ < e for some y G A}, and ||-|| is the Euclidiean norm. Also, 



(3(P,Q) = sup 



j fd(P-Q) 



BL 



< 1 , (5) 



where = 1 1 - 1 1 ^ + ||*|loo' anc ^ II 'Hi * s ^ e usua l Lipschitz seminorm (see [2J p. 394] for 

further details). The fact that we focus on the Prokhorov and the Fortet-Mounier metric 
is due to the following fact, proved in [U Th. 11.7.1]. For any two sequences {Pi} , {Qi} C 
P (M fc ), the following three conditions (A)-(C) are equivalent: (A) lim^ +00 j3 (Pi, Qi) = 
0; (B) lim^ +00 p (Pi, Qi) = 0; (C) on some auxiliary probability space (£2* , J 7 * ,F*), there 
exist sequences of random vectors {N* (I) : I > 1} and {I* (/):/> 1} such that 

C (I* (I)) = Pi and C (N* (/)) = Q { for every I, and (I) - N* (Z)|| -» 0, a.s.-P*, (6) 

where £ (•) indicates the law of a given random vector, and \\-\\ is the Euclidean norm. 
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3 Main results 

Fix integers k > 1 and d\, ...,dk > 1, and consider a sequence of &:-dimensional random vectors 
of the type 



/ > i 



(7) 



where, for each / > 1 and every j = 1, k, /, is an element of Sj &dj . We will suppose the 
following: 



There exists rj > such that 



0) 



> rj, for every j = 1, k and every I > 1. 



• For every j = 1, k, the sequence 



E 



U) 



dj\ 



, l>l, 



(8) 



is bounded. 



Note that the integers d\, dk do not depend on I. For every I > 1, we denote by N (/) = 
a centered A;-dimensional Gaussian vector with the same covariance matrix as 



I (/), that is, 



E 



E 



1 d i [fi 



p(0 



(9) 



for every 1 < i,j < k. For every A = (Ai,...,Afe) € we also use the compact notation: 

<a,i(0>* = V* and (A,N(0) fc = EL^f. 



SI 



The next result is one of the main contributions of this paper. Its proof is deferred to Section 



Theorem 1 Let the above notation and assumptions prevail, and suppose that, for every j 



1, 



, k, the following asymptotic condition holds: for every p = 1, cL- 



1, 



0, as Z — > +oo. 



fl 02(d J - P ) 

Then, as I — > +oo and /or every compact set M C 

sup |E [exp (i (A, I (!)>*)] - E [exp (i (A, N (/)),) 
agm 



0. 



(10) 



(11) 



We now state two crucial consequences of Theorem [TJ The first one (Proposition [2]) provides 
a formal meaning to the intuitive fact that, since (jlip holds and since the variances of 1(1) do 
not explode, the laws of 1(1) and N (I) are "asymptotically close". The second one (Theorem 
[3]) combines Theorem [1] and Proposition [2] to obtain an exhaustive generalization "without 
covariance conditions" of Theorem (see the Introduction). Note that in the statement of 
Theorem [3] also appear Malliavin operators, so that our results are a genuine extension of the 
main findings by Nualart and Ortiz-Latorre in |13j . We stress that multiple stochastic integrals 
of the type I d (/), d > 1 and / G Sj ed , are always such that I d (/) G D 1,2 . 



5 



Proposition 2 Let the assumptions of Theorem [7] prevail (in particular, U0\) holds), and de- 
note by C (I (Z)) and £(N(Z)) ; respectively, the law of 1(1) and N(Z), / > 1. Then, the two 
collections {£(N(/)) : Z > 1} and {£(I(Z)) : Z > 1} are tight. Moreover, if 7 (•, •) metrizes the 
weak convergence on P (R k ) , then 



lim 7 (£(I(Z)),£(N(Z))) = 0. 

I— >+oo 



(12) 



Proof. The fact that {C (N (/)) : Z > 1} and {£ (I (/)) : / > 1} are tight is a consequence of 



the boundedness of the sequence 



and of the relation E[I^. (f l 



0')\2l 



E[(N, 



. The rest of 



the proof is standard, and is provided for the sake of completeness. We shall prove (|12p by contra- 
diction. Suppose there exist e > and a subsequence {/„} such that 7 (C (I (/„)) , C (N (l n ))) > £ 
for every n. Tightness implies that {l n } must contain a subsequence {l n '} such that C (I (l n ')) 
and £(N(l n >)) are both weakly convergent. Since (fTT|) holds, we deduce that C (I (l n >)) and 
C (N (l n ')) must necessarily converge to the same weak limit, say Q. The fact that 7 metrizes 
the weak convergence implies finally that 



7 (C (I (l n ,)) , C (N (| n /))) < 7 (£ (I (/„,)) , Q) + 7 (£ (N (/„>)) , Q) 



0. 



(13) 



thus contradicting the former assumptions on {/ n } (note that the inequality in (fT3j) is just the 
triangle inequality). This shows that (|12p must necessarily take place. ■ 

Remarks, (i) A result analogous to the arguments used in the proof of Corollary [2] is stated 
in [U Exercise 3, p. 419]. Note also that, without tightness, a condition such as (lllj) does 
not allow to deduce the asymptotic relation (|12p . See for instance [U Proposition 11.7.6] for a 
counterexample involving the Prokhorov metric on P(R). 

(ii) Since (|12p holds in particular when 7 is equal to the Prokhorov metric or the Fortet- 
Mounier metric (as defined in (j4]) and ([5])), Proposition [2] implies that, on some auxiliary 
probability space (0* , T* ,P*), there exist sequences of random vectors {N* (I) : I > 1} and 
{I* (/):/> 1} such that 



I* (/) l = I (Z) and N* (0 l = N (/) for every /, and ||I* (I) - N 



where 



-» 0, a.s. 

stands for the Euclidean norm (see ©, as well as [31 Theorem 11.7.1]). 



(14) 



Theorem 3 Suppose that the sequence 1(1), I > 1, verifies the assumptions of this section (in 
particular, for every j = l,...,k, the sequence of variances appearing in |2j) is bounded). Then, 
the following conditions are equivalent. 



1. As I — > +00, relation U0\) is satisfied for every j = 1, k and every p = 1, dj — 1; 
2. 



lim p(£(I(Z)),£(N(Z))) 



lim /?(£(I(Z)),£(N(Z))) = 



(15) 



where p and (3 are, respectively, the Prokhorov metric and the Fortet-Mounier metric, as 
defined in and f5|); 



3. As I — > +00, /or every j 
4" 



1, Zc, 



E 



^ /1 



3E 



r(j) 



E 



3(d,!) 2 



r(j) 
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4- For every j = 1, k, 



^/?(r(/,,, (v; 



a) 



0, (16) 



where p and (3 are the Prokhorov and Fortet- Mounter metric on R; 
5. For every j = 1, k, 



D 



Id, If, 



r(?) 



dj(d 3l 



sf> d i 



0, in L 2 (X) , 



(17) 



as I — > +oo, where D is the Malliavin derivative operator defined in Section^ 



Proof. The implication 1. ==> 2., is a consequence of Theorem Q] and Proposition [2j 
Now suppose (|15p is in order. Then, according to [J[ Theorem 11.7.1], on a probability space 
(0*,JF*,P*), there exist sequences of random vectors N* (I) = (N*' {1 \ N*' ij) ), I > 1, and 



I* (/) 



r*.(l) r*>(*0 



I 



I > 1, such that (|14p takes place. Now 



3E 



T di f, 



r(j) 



1 2 



3E 



(.7) 



E 



for every j = 1, k, so that 

i" 



E 



r(j) 



'I 



3E 



E* 



r*.0') 



E* 



2V, 



/v, 



0. 



(18) 



The convergence to zero in (fTHj) is a consequence of the boundedness of the sequence (jSJ), implying 
that the family = 



r*j(?) 



/V, 



I > 1, is uniformly integrable. To see why {^4^} is 

uniformly integrable, one can use the fact that, since each has the same law as an element 
of the djth chaos of X and each N?'^' is Gaussian, then (see e.g. Ch. VI]) for every p > 2 
there exists a universal positive constant C p j (independent of I) such that 

i/p 



E[|^*| p ] 1/p 



= E* 
< E* 



r*.C?) 



*>(j) 



dp' 



4/4p 



+ E 



2V, 



*>(i) 



ip' 



4/4p 



< (7 PJ -E* 



r*>C?) 



N, 



*>0') 



2C pJ x (djl) 



f, 



U) 



where Mj = sup^ (djl) < +oo, due to (JSj) - This proves that 2. =>■ 3.. The implication 

3. ==> 1. can be deduced from the formula (proved in |14[ p. 183]) 

n 2 



E 



r(») 



3E 



2 dj /, 



E 



r(») 



2( dj -p) 



+ 



2{d J -p) 
dj-p 



Sf> 2 { d 3 



-)}' 
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The equivalence 1. -4=>- 4. is an immediate consequence of the previous discussion. 

To conclude the proof, we shall now show the double implication 1. <J=^ 5.. To do this, we 
first observe that, by performing the same caclulations as in |13[ Proof of Lemma 2] (which are 
based on an application of the multiplication formulae for multiple integrals, see [HI Proposition 
1.1.3]), one obtains that 



D 



hAft 



r(j) 



djidMfi 



+4 ]>>-!)! 

P =i 



h {d] - P ) (fi 3) ® P fl j) 



®p ft 



U) 



> 



2(d j - P ) 



the last relation implies immediately that 



Since 

1. 5.. To prove the opposite implication, first observe that, due to the boundedness of 
and the Cauchy-Schwarz inequality, there exists a finite constant M (independent of j and /) 
such that 



f! 3) ® P fP 



i{d j -p) 



< 



(?) 



< M. 



This implies that, for every sequence {/ n } ; there exists a subsequence {l n '} such that the se- 



quences 



fulfil; 



Ai) 



are convergent for every j = 1, k and every 



(?) 



$i j 



p = 1, —,dj — 1 (recall that, by assumption, there exists a constant r] > 0, such that f[ ; 
> rj, for every j and I). We shall now prove that, whenever (|17p is verified, then necessarily 
Jl ,®pJl , ®2(d,- P ) ~^ 0- Indeed, Theorem 4 in [T3] implies that, if (|TT|) takes place and 



Sf> d 3 



— » c > 0, then necessarily 



iV(0,c) 



(19) 



where iV (0, c) stands for a centered Gaussian random variable with variance c. But Theorem 1 



in [U] implies that, if (USD is verified, then f^ <g> p /, u ' 
This shows that 5. =>- 1.. ■ 



(?) 



0, thus proving our claim. 



The next result says that, under the additional assumption that the variances of the elements 
of I (I) converge to one, the asymptotic approximation (fT5|) is equivalent to the fact that each 
component of I (I) verifies a CLT. The proof is elementary, and therefore omitted. 

Corollary 4 Fix k > 2, and suppose that the sequence 1(1), I > 1, is such that, for every 
j = l,...,k, the sequence of variances appearing in |3|) converges to 1, as I — > +oo. Then, each 
one of Conditions 1.-5. in the statement of Theorem^ is equivalent to the following: for every 
j = 1, k, 

Id, iff) ^ #(0,1), (20) 
where N (0, 1) is a centered Gaussian random variable with unitary variance. 
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Remark. The results of this section can be suitably extended to deal with the Gaussian 
approximations of random vectors of the type (F^ (X) , (X)), where F^'(X), j = 
l,...,k, is a general square integrable functional of the isonormal process X, not necessarily 
having the form of a multiple integral. See [10\ Th. 6] for a statement containing an extension 
of this type. 



4 Proof of Theorem [T] 

We provide the proof in the case where 

Sj = L 2 ([0, 1] ,B ([0, 1]) , dx) = L 2 ([0, 1]), (21) 

where dx stands for Lebesgue measure. The extension to a general 9) is obtained by using the 
same arguments outlined in |14^ Section 2.2]. If fj is as in (I2ip . then for every d>2 one has that 
Sj Qd = L 2 ([0,l] d ), where the symbol L 2 ([0, l] d ) indicates the class of symmetric, real-valued 
and square-integrable functions (with respect to the Lebesgue measure) on [0,1] d . Also, the 
isonormal process X coincides with the Gaussian space generated by the standard Brownian 
motion 

t»W t ±X (l [0)t] ) , t E [0, 1] . 

This implies in particular that, for every d > 2, the Wiener-Ito integral I d (/), / E L 2 S ([0, l} d 
can be rewritten in terms of an iterated stochastic integral with respect to W, that is: 

h(f) = d\ I f 1 ■■■ f(t 1 ,...,t d )dW td ---dW t2 dW tl . (22) 

Jo Jo Jo 

We also have that h (/) = £ / (s) dW s for every / E L 2 ([0, l] 1 ) = 1? ([0, 1]). Note that the 
RHS of (|22p is just an iterated adapted stochastic integral of the Ito type. Finally, for every 
/ E L 2 ([0, l} d ), every g E -^([0, l] d ) and every p = 0, d A d', we observe that the contraction 
/ ® p g is the (not necessarily symmetric) element of L 2 ([0, \\ d+d ~ 2p ) given by: 

f® P 9 (yi, ...,yd+d'-2 P ) = / f{yi,->yd-p,ai,...,ap)x (23) 

^[o,i] p 

x 9 {Vd-p+h ■■■,yd+d'-2 P ,ai, ...,a p ) dax...da p . 

In the framework of pip , the proof of Theorem Q] relies on some computations contained 
in |15j . as well as on an appropriate use of the Burkholder-Davis-Gundy inequalities (see for 
instance |16l Ch. IV §4]). Fix A = (Ai, —,\k) G an d consider the random variable 

k 



^ PL P u d ■ —1 

(X,l(l)) k = ^Xjdjl •••/ 3 fi j) {u 1 ,...,u dj )dW Udj ---dW Ul 

Jo Jo 

A//,!-/,), (if) = f [EVilJS-i {fi j) («»■)) 



k 

E- 



dVF„. 
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where, for every d > 1, every t G [0,1] and every f € L% ^[0, l] d j, we define J l d {f) = 

Id ^/lj ot jd^ jd\ (for every c 6 1, we also use the conventional notation J\ (c) = c). We start 
by recalling some preliminary results involving Brownian martingales. Start by setting, for 
every u G [0,1], ^ («) = £*=! W*-l (f^ ( 



application 



3=1 J ° 



1| u j(j r i ), and observe that the random 



<t>XJ («) t G [0, 1] , 



defines a (continuous) square-integrable martingale started from zero, with respect to the canon- 
ical filtration of W, noted {j-f : t G [0, 1]}. The quadratic variation of this martingale is classi- 
cally given by 1 1— > Jq 0aj (^) 2 an d a standard application of the Dambis, Dubins and Schwarz 
Theorem (see |16t Ch. V §1]) yields that, for every I > 1, there exists a standard Brownian 

motion (initialized at zero) such that 



(A, 1(0), 



Jo 



cf>x,i (u) dW u 



W ( \' l) 2 

Jo du 



Note that, in general, the definition of W^'^ strongly depends on A and /, and that W^' 1 ^ is not 
a J^-Brownian motion. However, the following relation links the two Brownian motions W^ X '^ 

and W: there exists a (continuous) filtration |g t (A '° : t > o} such that (i) W t {X,l) is a Q[ X ^- 
Brownian motion, and (ii) for every fixed s G [0, 1] the positive random variable J Q S <f>\j (u) 2 du 



is a C/^'^-stopping time. Now define the positive constant (which is trivially a (/^''"'-stopping 
time) 



(x,l) 



q(X,l) 



E(0 A)i (u) )du, 



and observe that the usual properties of complex exponentials and a standard application of the 
Burkholder-Davis-Gundy inequality (in the version stated in [161 Corollary 4.2, Ch. IV ]) yield 
the following estimates: 



E[exp (i(A,I(0) fe )]-E[exp (iW ( 



(A,0 

q(X,l) 



< 



exp ( iwW . . a . 



E 



exp ( iW 



r(A,0 



ff(A,i) 



/o 4>x,l(u) du 



W 



(X,l) 



q(X,l) 



< E 



< CI 



/0 </!j,lW d « 



(A,I) 



q(X,l) 



<Px,i (u) 2 du -q(X,l) 



(24) 



where C is some universal constant independent of A and To see how to obtain the inequality 
(f24l) . introduce first the shorthand notation T (A, I) = (f)\j (u) 2 du (recall that T(X,l) is a 



? (A,0 



stopping time), and then write 



J,, 1 <t>\,i( u ) 2 du 



w 



(X,l) 



q(X,l) 



T(X,l)Vq(X,l) 



T(X,l)Aq(X,l) 



dW^ l) 



T(X,l)Vq(X,l) 



H (u) dW^ x ' l) 
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where H (u) is the (^'^-predictable process given by H (u) = 1 {u > T (A, I) A q (A, /)}, so that 
rT(A,0Vg(A,0 



-ff (ti) du 



|T(A,/)Ag(A,/)-T(A,/)Vg(A,0| 



</>A,K n ) du-q(X,l) 



= \T(X,l)-q(X,l)\ 

In particular, relation (|24p yields that the proof of Theorem [1] is concluded, once the following 
two facts 
sequence 



two facts are proved: (A) = (A, N(Z)) fc , for every A G M fc and every I > 1; (B) the 



E 



0AJ (u) 2 du - q (X, I) 



I > 1, 



converges to zero, uniformly in A, on every compact set of the type M = [— T, T] k , where 
T £ (0, +oo). The proof of (A) is immediate: indeed, W^ X ' 1 ^ is a standard Brownian motion 
and, by using the isometric properties of stochastic integrals and the fact that the covariance 
structures of N (I) and I (I) coincide, 













Qf * ^ (u) dW u \ 


= E 




= E 


jA,N(Z)>J 













(A,/) = / E(<j> Xjl (uY)du = E 



To prove (B), use a standard version of the multiplication formula between multiple stochastic 
integrals (see for instance [12[ Proposition 1.5.1]) 



l * 



i / k 

l Ew«-i(/i 



i,<=l 



/ j) («r)l [0)U] <i f -i 



dw, 



'(A,0 + 



E***jfg(V)(V) 



(25) 



where the index D (i,j) is defined as 



D(i,j) 



di - 2 



if di = di 



min (di , dj) — 1 if dj 7^ dj . 
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Formula (|25p implies that, for every A G [-T,T] k (T > 0), 



E 



</> A j( n ) du-q(X,l) 



k D(i,j) . . . 



(26) 



xE 



idi+dj-2-2p ((/, W («, •) Ip.u]'*!- 1 ) ®P (fl J> ( u i •) 1 p,u] ,, i- 1 ) 



(i) 



du 



(note that the RHS of (i26|) does not depend on A) . Finally, a direct application of the calculations 
contained in [151 P- 253-255] yields that, for every i,j = 1, k and every p = 0, D 



E 



0, (27) 



as Z — > +oo. This concludes the proof of Theorem [TJ 



Remark. By inspection of the calculations contained in ^5 ( p. 253-255], it is easily seen 
that, to deduce ([27]) from (fTUj) . it is necessary that the sequence of variances (|8|) is bounded. 



5 Concluding remarks on applications 

Theorem [1] and Theorem [3] are used in [10] to deduce high-frequency asymptotic results for 
subordinated spherical random fields. This study is strongly motivated by the probabilistic 
modelling and statistical analysis of the Cosmic Microwave Background radiation (see [7], [8], 
[9] and |10j for a detailed discussion of these applications). In what follows, we provide a brief 
presentation of some of the results obtained in |10j . 

Let S 2 = {i £ I 3 : ||x|| = l} be the unit sphere, and let T = {T (x) : x G S 2 } be a real- 
valued (centered) Gaussian field which is also isotropic, in the sense that T (x) = T (Tlx) (in 
the sense of stochastic processes) for every rotation 1Z £ SO (3). The following facts are well 
known: 

(1) The trajectories of T admit the harmonic expansion T (x) = ][^ Ylm=-i a imXim where 

\Xlm '■ I > 0, m = —/,...,/} is the class of spherical harmonics (defined e.g. in [17, Ch. 

5]); 

(2) The complex-valued array of harmonic coefficients {ai m : I > 0, I > 0, m = is 

composed of centered Gaussian random variables such that the variances E|aj m | 2 = Ci 
depend exclusively on I (see for instance [I]); 

(3) The law of T is completely determined by the power spectrum {C/ : I > 0} defined at the 

previous point. 
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Now fix q > 2, and consider the subordinated field 

T^(x)±H q (T(x)), xeS 2 , 

where H q is the gth Hermite polynomial. Plainly, the field TW is isotropic and admits the 
harmonic expansion 

oo I oo 

Tiq) w = EE a '-^ ( x ) - E T / 9) (*) > 

1=0 m=-l 1=0 

where ai m - q = / §2 T (<7) (z) Yi m (z)cte. For every Z > 0, the field t/ 9) = Yl m =-i a im;qYi m is real- 
valued and isotropic, and it is called the Zth frequency component of (see [7J or [10 for a 
physical interpretation of frequency components). In [10J, the following problem is studied. 

Problem A. Fix q > 2. Find conditions on the power spectrum {Ci : I > 0} to have that 
the finite dimensional distributions (f.d.d.'s) of the normalized frequency field 

Var (T l {q) (x) 



are 'asymptotically close to Gaussian', as I — * +oo. 

The main difficulty when dealing with Problem A is that (due to isotropy) one has always 
that 

E[T[ 9) (x)Ti 9) (y)] =P l (co S (x,y)), (28) 

where Pi is the Zth Legendre polynomial, and (x, y) is the angle between x and y. Indeed, since 
in general the quantity Pi (cos{x,y)) does not converge (as I — ► +oo), one cannot prove that 

the f.d.d.'s of t[ 9 ^ converge to those of a Gaussian field (even if (x) converges in law to 
a Gaussian random variable for every fixed x). However, as an application of Theorem Q] and 
Proposition [21 one can prove the following approximation result. 

Proposition 5 Under the above notation and assumptions, suppose that, for any fixed x E § 2 , 

T^ix) L -T AT (0,1). (29) 

i^+oo 



Then, for any k > 1, any xi,...,Xk G 8 2 and any 7 metrizing the weak convergence on P 

1 (c(T { l q) (x 1 ),...,T i l g) (x k )),N(l)) ^0, (30) 

where, for every I, N(i) = (n^\...,N^ is a centered real-valued Gaussian vector such that 

E {jV, W JV, y) } = ^ (cos (xi, Xj)) . 
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Proof. Since (x) is a linear functional involving uniquely Hermite polynomials of order 
q (written on the Gaussian field T) one deduces that there exists a real Hilbert space f) such 
that (in the sense of stochastic processes) 

(X) I q (f (gM ) , 

where the class of symmetric kernels 

{/(,,!,*) : I > 0, x e S 2 } 

is a subset of fj 013 , and I q (f(q,i,x)) stands for the qth Wiener-Ito integral of f( q ,i, x ) with respect 
to an isonormal Gaussian process over Sj, as defined in Section [2j Since the variances of the 
components of the vector (T[ 9 (x\) , t\ 9 ^ (xk)) are all equal to 1 by construction, we can apply 
Theorem [3] and Proposition [2 Indeed, by Theorem [3] we know that (I29jl implies that, for every 
p = 1, q — 1 and every j = 1, k, 

f( q ,l, x ,)® P f( q ,l, Xj )^0m 

Finally, Proposition [2] and (|28p imply immediately the desired conclusion. ■ 

The derivation of sufficient conditions to have (|29p is the main object of [10] . In particular, 
it is proved that sufficient (and sometimes also necessary) conditions for (|29p can be neatly 
expressed in terms of the so-called Clebsch-Gordan coefficients (see again [H]), that are elements 
of unitary matrices connecting reducible representations of SO (3). 
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